DISTRIBUTION OF RESONANCE WIDTHS AND DYNAMICS OF CONTINUUM 
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We analyze the statistics of resonance widths in a many-body Fermi system with open decay 
channels. Depending on the strength of continuum coupling, such a system reveals growing devi- 
ations from the standard chi-square (Porter-Thomas) width distribution. The deviations emerge 
from the process of increasing interaction of intrinsic states through common decay channels; in 
the limit of perfect coupling this process leads to the super-radiance phase transition. The width 
distribution depends also on the intrinsic dynamics (chaotic vs regular). The results presented here 
are important for understanding the recent experimental data concerning the width distribution for 
neutron resonances in nuclei. 
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Open and marginally stable quantum systems are of 
great current interest in relation to numerous applica- 
tions in nuclear physics of exotic nuclei, chemical reac- 
tions, condensed matter, astrophysics and quantum in- 
formatics. The general problem can be formulated as 
that of the signal transmission through a complicated 
quantum system. The complexity of theoretical descrip- 
tion of such processes is due to the necessity of a consis- 
tent unified theory that would cover intrinsic structure, 
especially for many-body systems, along with cross sec- 
tions of various reactions. 

One of the best and historically advanced examples of 
the manifestation of the interplay between intrinsic dy- 
namics and decay channels is given by low-energy neu- 
tron resonances in complex nuclei [l[ . The series of these 
well pronounced separated resonances were studied long 
ago 043 and later gave rise to the "Nuclear Data En- 
semble" d, @. Interpreting these resonances as quasi- 
stationary levels of the compound nucleus formed after 
the neutron capture, agreement was found with predic- 
tions of the Gaussian Orthogonal Ensemble (GOE) of 
random matrices. With exceedingly complicated wave 
functions of compound states, the statistical distribution 
of their components is close to Gaussian. The neutron 
decay implements the analysis of a specific component re- 
lated to the channel "neutron in continuum plus a target 
nucleus in its ground state" . The neutron width is pro- 
portional to the squared amplitude of this component, 
and the width distribution then is xt with v = 1 as 
appropriate for one channel [Porter- Thomas distribution 
(PTD)]. 

The recent experiments with improved accuracy Q 
give evidence of significant deviations from the PTD so 
that the attempts to still use the xl distribution for the 
fit invariably require v < 1. A non-pure set of resonances, 



for example, a sequence of mainly s-resonances contam- 
inated by p-wave states, would shift the distribution to 
a higher number of degrees of freedom, v > 1. The new 
result was interpreted as a consequence of an unknown 
non-statistical mechanism or just a breakdown of nuclear 
theory as was claimed in the related article in "Nature" 
with a title "Nuclear theory nudged" 0. 

The goal of this letter is to point out that a correct 
description of unstable quantum states in a complicated 
many-body system naturally leads to deviations from the 
GOE and PTD, of the same type as observed in Q. The 
random matrix theory was formulated for local statis- 
tics in a closed quantum system with a discrete spec- 
trum governed by a very complicated Hermitian Hamil- 
tonian. As such, its predictions were repeatedly checked, 
qualitatively and quantitatively, in systems like quantum 
billiards Q and their experimental embodiment in mi- 
crowave cavities \WL_ Tl| , and in shell-model calculations 
for complex atoms [12j and nuclei (l3j . 

However, the presence of open decay channels and 
therefore the finite lifetime of intrinsic states unavoid- 
ably lead to new phenomena outside of the GOE frame- 
work [3, [TH as it was clearly demonstrated by the first 
numerical simulations |16h18| . Two interrelated effects 
follow from the fact that we deal with unstable rather 
than with strictly stationary states: the level repulsion 
disappears at the spacings comparable to the level widths 
and the growing widths undergo the redistribution with 
the trend to collectivization and eventually formation of 
super-radiant (short-lived) states along with the narrow 
(trapped) states [H, HU . The new dynamics modify the 
GOE predictions as well as certain features of Ericson 
fluctuations [20( in the regime of overlapping resonances 
[2l| . The occurrence of a super-radiant transition has 
been also demonstrated outside the random matrix the- 
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ory framework f22| . 

A quantum many-body system coupled to open decay 
channels can be rigorously described by an effective non- 
Hermitian Hamiltonian T-L [23| , 

n = H- % -w. (i) 

Here H is the Hcrmitian Hamiltonian of the closed sys- 
tem that in general includes virtual (off-shell) coupling 
to the continuum, while the anti-Hermitian (on-shell) 
part (—i/2)W is constructed in terms of the amplitudes 
A = {Af} coupling intrinsic states \i) to the open chan- 
nels c, 

M 

W = AA T => Wij = A t A r ( 2 ) 

c(opcn) 

We consider a time-reversal invariant system, when the 
(in general depending on running energy) amplitudes A1 
can be taken real. The factorized form of W is an impor- 
tant property that follows from unitarity of the scattering 
matrix [241 ]. The complex eigenvalues, £ — E— (i/2)T, of 
the Hamiltonian (fTJ) coincide with the poles of the scat- 
tering matrix and determine the positions and the widths 
of the resonances in cross sections of various reactions. 

The factorized matrix W has M non-zero eigenval- 
ues, this number being equal to that of open channels. 
The matrix H has dimension N that in the nuclear case 
should include a large number of shell-model many-body 
states important for the dynamics in the energy range 
under consideration; in the region of neutron resonances, 
N ~ 10 5 ^ 6 > M. With the trace of W equal to 77, the 
important parameter defining the dynamics is the ratio 
of typical "bare" widths rj/N of individual states to the 
energy spacings D. At small value of this parameter, an 
open channel serves as an analyzer that singles out a spe- 
cific component of the exceedingly complicated intrinsic 
wave function. Characteristically, the resonance widths 
in such a system obey the PTD. With widths increas- 
ing, the system moves in the direction of the regime of 
overlapping resonances. 

Let us consider the single-channel case, M = 1, having 
in mind the s-wave elastic neutron resonances. With a 
high level density of intrinsic states at relevant energy, 
their local spectral statistics is close to the predictions 
of the GOE. Then, essentially owing to the central limit 
theorem, the individual components of a typical intrinsic 
state are Gaussian distributed uncorrelated quantities, 
and the neutron widths, being proportional to the abso- 
lute magnitudes of those components, display the PTD. 
However, the correct description of the dynamics with 
the continuum coupling shows the limited character of 
this prediction. The imaginary part ^ works similar 
to the collective multipole forces and creates the interac- 
tion between intrinsic states through continuum. When 
the coupling is weak, k = tj/ND < 1, we indeed ex- 
pect to see well isolated resonances with the PTD of the 
widths. With growing continuum coupling (increase of 



energy from the threshold), the deviations become more 
and more pronounced. At k ~ 1, a kind of a phase tran- 
sition occurs with the sharp redistribution of widths and 
the segregation of a super-radiant state accumulating the 
lion's share of the whole summed width, an analog of a 
giant resonance along the imaginary energy axis. The 
formal mechanism is clear from the factorized structure 
of W that, for M = 1, has only one non-zero eigenvalue 
equal to the trace of W. Being similar to super-radiance 
in optics [25l | , this mechanism works essentially indepen- 
dently of the regular or chaotic nature of intrinsic dy- 
namics. 

In the case of the GOE-type dynamics of the closed 
system and M = 1, the distribution of complex eigenval- 
ues can be found analytically. The Ginibre ensemble [26[ 
of complex Gaussian matrices is not applicable since it 
considers the imaginary parts of the eigenvalues spread 
over the complex plane while physical widths are posi- 
tive. The exact result, see derivation in [15j . is given by 
the Ullah distribution (27|, 

V({E n ;T n }) = C N H ^ ~ J] 4= cxp(-AF). 

(3) 

Here Cjv is a normalization constant; the pre-exponential 
factor describes the correlations which are reduced to the 
usual GOE level repulsion for stable states but for com- 
plex energies £ n contain the interactions with their "elec- 
trostatic images" £ * . Along with the Porter-Thomas fac- 
tor r^ 1 / 2 , this guarantees that the widths are positive. 
The "equilibrium" distribution of complex energies is de- 
termined by their "free energy" , 

n x ' m<n 

that includes the interaction between the widths, the 
last term in Eq. (j4j. The mean level spacing in the 
closed system is defined as D = 2a/N, where 2a is the 
spectral interval of real energy. The Gaussian ensemble 
of the decay amplitudes is defined by the mean values 
A n = 0, A n A m = (r}/N)6mn- Here, a regular evolution 
of the widths as a function of energy of the resonance is 
excluded as it is usually done with the rescaling to the re- 
duced widths; we do not discuss here the way of practical 
rescaling that may depend on the specific nucleus. Our 
main purpose here is to show that systematic deviations 
from the PTD occur even for the set of reduced widths, 
and the effects are caused by the interaction (0}. 

For very small widths, k <C 1, the width interaction 
is negligible, the first product in © reduces to the stan- 
dard level repulsion, and the distribution Q is factorized 
into the product of the GOE distribution of real energies 
and the PTD for the widths. While the usual Hermi- 
tian perturbation mixes the intrinsic wave functions and 
therefore makes their widths close to each other (level 
repulsion and width attraction), the anti-Hermitian in- 
teraction through the continuum relaxes the level repul- 
sion but leads to the collectivization through the common 
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decay channel and width repulsion (28|. In the case of 
k ~ 1, the width repulsion becomes critical, and the most 
probable configuration is the one where this repulsion is 
small because the total width is going to concentrate in 
a single super-radiant state [HI, [l9| . With the further- 
increase of k 3> 1 (but for a fixed number of open chan- 
nels), the broad state becomes a smooth envelope, and 
we return to the set of (TV — 1) narrow resonances. Below 
we show the typical evolution of the width distribution 
studied in a large-scale numerical simulation. 

In the first type of simulations, we considered H as a 
member of the GOE [3(| Hl| where the matrix elements 
are Gaussian random variables, (Hfj) — (l+5ij)/N. This 
corresponds to the limiting case of fully chaotic intrin- 
sic dynamics. In parallel, we also modeled H by the 
two-body random ensemble (TBRE) for n fermions dis- 
tributed over m single-particle states; the total number 
of many-body states is N — m\/[n\{m — n)!]; in our sim- 
ulations n = 6, m = 12, N = 924. The TBRE is 
modeled by the Hamiltonian H = Hq + V, where the 
mean field part Ho is defined by single-particle energies, 
ej, with a Poissonian distribution of spacings and the 
mean level spacing do- The interaction V [29[ is fixed by 
the variance of the two-body random matrix elements, 
(Vi2-3 ,4) = v o- While at v — we have a Poisso- 
nian spacing distribution P(s) of many-body states, for 
do = (infinitely strong interaction, vo/do — >• oo), P(s) 
is close to the Wigner-Dyson distribution typical for a 
chaotic system. The critical interaction for the onset of 
strong chaos is given [29j by v cr /do ps 2(m—n)/N s , where 
N s = n(m — n)[l + (n — l)(m — n— l)/4] is the number 
of directly coupled many-body states in any row of the 
matrix Hij. In our model, v cr /do ~ 1/20. 

We computed the complex eigenvalues of the effective 
Hamiltonian, Eq. ([!]), for 10 3 random realizations, select- 
ing the states at the center of the energy band where 
their density is almost constant. 

The distribution of widths, normalized to their average 
value, was obtained for different strengths of continuum 
coupling. The numerically obtained distributions were 
fitted, using a standard x 2 test, with a xt distribution 
for v degrees of freedom, following a common practice in 
nuclear data analysis. As a measure of quality of the fit 
we used the criterion x 2 ~ 1> where x 2 is the reduced 
chi-square value [32| . Fig. [T] shows the normalized width 
distribution for the GOE case with M = 1, left pan- 
els, and M = 2, right panels, for different coupling to 
the continuum, n. The standard distributions, PTD for 
M = 1, and x1=2 f° r M = 2, are valid only when the 
coupling to the continuum is very weak, while strong de- 
viations from xt=M appear, both for large and small T, 
as we increase the coupling. In Fig. [T] we also show the 
best fit possible for a x 2 distribution; always the corre- 
sponding value of v is < M, moreover the quality of the 
fit decreases as k increases, see discussion below. In the 
lowest panels in Fig. [1] the distribution of the widths is 
shown at the super-radiance transition, k = 1. The tail 
here is described by a power law, see discussion in Ref. 



[33|, meaning that no xl distributions would be a good 
fit. 



M=l M=2 




0.1 1 10 0.1 1 10 

r/<r> r/<r> 



FIG. 1: (Color online) Normalized distribution of the widths 
for one open channel, M = 1, left panels, and for M = 2, 
right panels, for the GOE intrinsic Hamiltonian and different 
continuum coupling strength, k. The numerical results are 
shown by histograms; the expected PTD for M — 1 and xt=2 
distribution for M = 2 are shown by a smooth curve in all 
panels. With full squares we show the best fit with a xt 
distribution. In the lowest panels, a straight line shows the 
power law P(T/{T}) oc (r/(r))~ 2 . Included are the widths 
for the resonances in the energy interval ±0.5 for 10 3 different 
ensemble realizations. 

The dependence of the best fitted value of v on the 
coupling strength is shown in Fig. [5] for M — 1, upper- 
panel, and for M = 2, lower panel. Along with the data 
for the GOE intrinsic Hamiltonian (circles), the data for 
the TBRE intrinsic dynamics are shown by crosses for 
the case vq = do/50 and by squares for vo — (no in- 
trinsic chaos). Regardless of intrinsic dynamics, the best 
value of v steadily decreases as n increases. It is clear 
that a family of x 2 - distributions is not appropriate to 
fit our numerical data, except for very weak continuum 
coupling. Indeed, the Xr criterion steeply increases with 
the coupling strength, see insets in Fig. [5J For weak in- 
ternal chaos, the departure from the xi=M distribution 
is stronger than for chaotic intrinsic dynamics, even at 
weak continuum coupling. The absence of intrinsic chaos 
and corresponding level repulsion implies a stronger sen- 
sitivity to continuum coupling. 

When analyzing empirical neutron s-wave resonances 
one of the main difficulties is the p-wave contamination. 
In order to analyze this problem we considered the case 
of two non-equivalent channels, with 10% of the states 
coupled to the additional channel by a smaller coupling 
strength. When another channel is included, the best fit 
value of v at weak coupling is always larger than one; 
again this value decreases when the continuum coupling 
is growing. As expected, we also observed the evolu- 
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FIG. 2: (Color online) The best fitted value of v vs the 
coupling strength to the continuum, k. Full circles refer to the 
case of GOE intrinsic Hamiltonian, crosses stand for the case 
of TBRE with vo = do/50, and squares for the case of TBRE 
with vo = (no intrinsic chaos); the number of channels is 
one (upper panel) and two (lower panel) . The insets show the 
reduced chi \r as a function of k. 

tion of the level spacing distribution P(s) along the real 
axis, with disappearance of repulsion at spacings compa- 
rable with the level widths [17|; for example, at n = 0.5, 
P(0) sa 0.2. Interesting physics related to energy-width 



correlations will be discussed elsewhere. 

To summarize, the normalized width distributions 
have been analyzed for a system with M — 1 or 2 
open channels as a function of the continuum coupling 
strength. As this coupling increases, the best fit value of 
v for a xt distribution decreases below M , in accordance 
to recent experimental findings Q. At the same time, the 
fit quality becomes poor showing that the standard PTD 
(and in general any xt distribution) is applicable only 
for extremely narrow resonances. The low-energy neu- 
tron resonances in a heavy nucleus correspond to the very 
beginning of the process of width collectivization. How- 
ever, already here the deviations from the (GOExPTD) 
factorized distribution are noticeable. These deviations 
are more pronounced for regular intrinsic dynamics than 
for chaotic intrinsic dynamics. Therefore, the interpreta- 
tion of the width as a strength of the pure neutron com- 
ponent in the compound wave function fails due to the 
coupling through continuum that has to be accounted for 
in a proper statistical description. The phenomenon un- 
der discussion is of general nature and it may influence 
all processes of signal transmission through a quantum 
system. 
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